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摘要 :秦岭 太白 山林 线 植被 因 海 拔 较 高 且 受 人 为 扰动 较 轻 , 对 气候 变化 的 响应 尤为 敏感 ,为 获取 过 去 气候 变化 信息 提供 了 可 靠 
代用 资源 。 然 而 ,结合 树木 年 代 学 方法 及 Arcgis 空间 插值 功能 进行 秦岭 林 线 气候 变化 重建 的 工作 宇 今 俏 处 科 空 白 。 利 用 采 上 
太白 山林 线 地 带 太白 红 杉 (Larix chinensis) 所 建立 的 树 轮 宽度 资料 ,与 提取 自 太白 山 保 护 区 气温 机 格 数据 中 的 采样 点 位 置 气象 
数据 进行 相关 分 析 。 结 果 表明 ,太白 红 杉 与 3 一 6 月 平均 气温 相关 性 最 显著 ,采用 线性 回归 建立 了 两 者 的 拟 合 模型 ,剔除 重建 方 
程 中 的 1997 .1998 年 之 后 ,方差 解释 量 达 57.2% (调整 自由 度 后 为 535.5% ) ;重建 气温 序列 显示 偏 冷 时 段 平 均 跨度 (16 年 ) 较 仿 
暧 时 段 平 均 跨度 (10.8 年 ) 长 , 偏 冷 时 段 有 :1870 一 1881 年 .1903 一 1918 年 和 1977 一 1996 年 ; 偏 暖 的 时 段 有 : 1882 一 1892 年 、 
1919 一 1929 年 和 1997 一 2013 年 ;在 1931 一 1978 年 这 一 时 期 ,气温 相对 稳定 ,1988 年 芍 后 升温 强烈 ;周期 分 析 显 示 近 144 年 以 来 
3 一 6 月 气温 存在 22 一 31 a,18 一 22 a 以 及 10 一 13 a 的 3 个 振荡 周期 ,可 能 与 大 尺度 气候 驱动 及 太阳 活动 存在 联系 。 以 上 结果 均 
得 到 历史 记录 以 及 周边 重建 结果 的 支持 。 
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Abstract: Since the vegetation of,the timberline area at Mount Taibai in the Qinling mountains is of high elevation and 
rarely disturbed by anthropogenic actiyities, it is sensitive to climate change, thereby providing a proxy resource of past 
climate change information\ Howeveérn, to date, there have been no investigations that have integrated dendrochronological 
methods with Arcgis spatial analysis techniques to reconstruct the past climate change at specific elevations and locations. 
Herein, we presented a réconstruction based on extracted meteorological data from spatial interpolation and Larix chinensis 
tree ring widths. Climaté growth response analysis revealed that the air temperature ,particularly from March to June, was the 
principlefactor Mthat, limited the radial growth of L. chinensis . By using dendrochronology techniques, we obtained a 
recolistruction of the March to June mean air temperature for the timberline of Mount Taibai over the period from 1870 to 
20138. The@ reconstruction captured 57.2% of the temperature variance( 55.5% after the degree of freedom is adjusted) after 
removing the tree—ring width values of 1997 and 1998 during the calibration periods. On the decadal scale, the pattern of 
cool=warm fluctuations exhibited frequent repetitive features., On average, the cold periods ( 16 years) always persist longer 
than the warm periods( 10.8 years).The years 1870 一 1881 ，1903 一 1918 and 1977 一 1996 were colder periods, whereas 
the years 1882 一 1892，1919 一 1929，1950 一 960 ，and 1997 一 2013 were relatively warmer. The temperature variations 
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indicated three quasi 一 periodic changes on scales of 22 一 31 a,18 一 22 a and 10 一 13 a, which suggested that the Mount 
Taibai region may be associated with a huge 一 scale climate drive and solar activity. Our reconstruction exhibited good 


consistency with local historical records and reconstruction results from similar climatic regions for the same intervals. 


Key Words: Qinling mountain; Timberline; Larix chinensis; Tree-ring width; Temperature interpolation; Climate 


reconstruction 


随 着 近年 来 全 球 气候 变化 对 生态 环境 的 影响 越 来 越 深 ,对 可 持续 发 展 构成 严重 威胁 ,已 经 被 世界 各 国 广 
泛 关注 员 。 为 了 应 对 未 来 气候 变化 的 不 确定 性 ,历史 时 期 气候 变化 规律 及 机 制 的 研究 价值 逐渐 显现 ,已 经 成 
为 PAGES 重点 研究 的 方向 之 一 站 。 然 而 , 因 地 势 环境 等 条 件 不 便 ,秦岭 地 区 高 山 气象 站 点 分 布 稀 路 ,气象 数 
据 相 对 短缺 ,大 多 数 台 站 资料 年 份 较 短 ,因此 利用 非常 规 的 代用 资料 进行 历史 气 俩 的 研究 近年 来 进展 十 分 迅 
速 ] 。 在 众多 代用 资料 中 ,树木 年 轮 与 其 他 代用 资料 相 比 具有 定年 准确 .连续 性 强 、 分 辩 率 穴 、 与 气候 要 素 相 
关 性 好 地域 分 布 广 旦 易于 采样 等 优点 ,从 而 得 到 了 广泛 的 应 用 "1 。 

作为 西北 地 区 力 至 全 国 重要 的 地 理 分 界线 ,在 过 去 的 20 年 内 ,秦岭 地 区 例 未 年 代 学 工作 发 展 十 分 迅 
速 ! 沾 。 刘 洪 斌 、 邵 雪 梅 “7 最初 重建 了 镇 安 地 区 1755 年 以 来 的 初春 温度 变化 ,并 发 现 镇 安 地 区 初春 温度 重建 
序列 具有 显著 的 准 50 a 及 准 2 一 3 a 周期 ,随后 结合 镇 安 鹰 嘴 崖 . 佛 坪 光 头 山 、 华 山 和 太白 山南 北 坡 年 表 , 重 建 
了 近 300 a 来 秦岭 地 区 的 初春 温度 变化 ;刘备 、 马 利 民 等 中 重建 有 过 去 251a 镇 安 地 区 3 一 4 月 平均 气温 序列 ， 
发 现 秦 岭 冷杉 与 该 地 区 早 涝 指数 存在 显著 负 相 关 关 系 , 可 视 作 坚 涝 灾情 的 预警 信号 ,而 刘备 、 刘 娜 ' 等 发 现 
了 秦岭 落叶 松 轮 宽 序列 与 1 一 7 月 平均 气温 存在 较 好 的 下 相关, 并 重建 了 该 地 区 1814 一 2003 年 1 一 7 月 的 平 
均 气 温 ; 戴 俊 虎 、 邵 雪 梅 "等 则 采用 太白 山高 山林 线 附近 的 太白 红 杉 树木 年 轮 资料 ,重建 了 太白 山地 区 近 300 
a 的 温暖 指数 ; 田 沁 花 中 等, 则 利用 伏 牛 山 油 松树 轮 宽 度 年 表 , 重 建 了 该 区 域 1874 年 以 来 5 一 7 月 平均 最 高 温 
度 , 发 现在 过 去 134 年 中 , 伏 牛 山区 5 一 7 月 平 的 最 高 气温 经 历 了 4 次 冷 期 和 5 次 暧 期 。 由 以 上 可 见 ,秦岭 地 
区 气候 重建 研究 工作 已 经 得 到 一 定 的 成 果 , 然 而 ,秦岭 以 东西 跨度 较 大 ,南北 地 形 复杂 著称 ,加 之 存在 气温 直 
减 率 的 客观 事实 ,重建 周边 气象 站 的 和 氢 候 变化 并 不 能 够 准确 代表 林 线 地 带 的 气候 变化 ,为 解决 这 一 问题 ,本 课 
题 组 尝试 借助 气象 站 资料 利用 Aregisi 软件 进行 空间 差 值 并 重新 提取 目标 位 置 海拔 的 逐年 气温 数据 以 重建 林 
线 地 带 的 气候 变化 过 程 ,研究 对 锡 太 白山 为 我 国 青藏 高 原 以 东 的 最 高 峰 。 

太白 山 作 为 全 球 和 区 域 气候 变化 的 敏感 区 域 '; ,虽然 气象 资料 短缺 ,但 太白 红 杉 作为 当地 唯一 的 可 成 纯 
林 的 树种 5 , 轮 宽 变 化 特征 明显 ' 提供 了 气候 变化 研究 的 理想 条 件 。 另 外 ,太白 山 为 西安 市 水 源 地 黑河 的 源 
头 汇 水 区 ,其 林 线 环境 的 变化 不 仅 对 周边 地 区 气候 产生 影响 ,还 关系 着 关中 ,西安 地 区 社会 和 经 济 发 展 ,重建 
其 林 线 地 区 气候 变化 具 明 确 指 示意 义 。 

本 文 以 酝 轮 样本 为 基础 ,以 气温 空间 持 值 技术 ' 趾 为 辅助 , 旨 在 准确 、 可 靠 得 重 建 近 144 年 太白 山林 线 地 
带 气温 变化 过 程 z 此 方法 在 秦岭 地 区 力 至 全 国 范围 内 仍 十 分 少见 ,不仅 能 丰富 西北 地 区 的 历史 气候 变化 资 
料 ,为 其 他 的 重建 资料 提供 对 比 数据 ,还 可 以 为 制定 太白 山林 线 生态 保护 决策 方案 提供 一 定 的 科学 依据 。 


DD 研究 区 概况 与 数据 处 理 


1.1 树 世 样本 的 采集 与 处 理 

研究 区 位 于 太白 山 最 高 点 拔 仙台 (3767 m) 附 近 的 太白 红 杉 森林 上 限 处 ,地 理 坐 标 为 (N 33°49' 一 34°10"， 
E107°19' 一 107°58') ,该 地 区 年 平均 气温 1.8 一 2.1% ,降雨 量 800 一 900 mm, 属 大 陆 性 季风 气候 区 ,冬季 多 风 少 
雨 ,夏季 高 湿 多 雨 。 太 白山 植被 带 垂 直 分 布 特点 明显 , 林 线 以 上 为 第 四 季 冰 期 形成 的 冰川 遗迹 和 高 山 草包 ,以 
下 接 巴 山 冷 杉 群 系 ,土壤 类 型 以 森林 草 旬 土 为 主 ' "1。 

本 课题 组 于 2013 年 10 月 沿 太 白山 东 侧 山 准 在 海拔 3235 m 的 太白 红 杉 纯 林 区 设立 了 采样 点 (图 1) , 树 
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木 平均 树 高 为 8.85 m, 平 均 胸 径 为 24.41 cm。 由 于 采样 
点 距离 上 板 寺 较 近 ,将 其 命名 为 SBS。 样 地 树 间 距 较 
小 ,树木 分 布 集 中 , 土 层 厚度 均 在 (20+10) em 范围 内 ， 
立地 坡度 均 在 (30+15)° 范 围 内 。 采 样 的 策略 依据 国际 


35°N 


树 轮 数据 库 的 标准 (ITRDB ) ,综合 了 敏感 性 原则 ,生态 34 | /7 J ee 
环境 原则 和 复 本 原则 等 ,挑选 健康 、 直 立 的 植株 进行 采 OY) 入 泰和 人 人 pm 人、 人 
样 ,使 用 生长 锥 在 胸 高 位 置 从 不 同方 向 取 A 、B 两 芯 作 和 山 及 A 

为 复 本 ,累计 采 得 50 个 原始 树 芯 样本 。 i 


树 轮 样 本 经 过 粘贴 、 了 晾晒 、 打 磨 等 一 系列 惯例 步 
又 预 处 理 后 ,随即 进行 初步 定年 工作 。 采 用 精度 为 0.01 
mm 的 LINTAB 轮 宽 仪 对 每 一 棵 树 世 每 一 年 的 树 轮 宽度 
进行 测量 。 完 成 测量 工作 后 使 用 COFECHA 程序 "中 对 
轮 宽 序列 进行 检验 和 交叉 定年 ,检验 过 程 中 ,对 受 外 界 
因素 导致 变异 的 树 芯 或 是 部 分 年 份 进行 筛选 .剔除 ,与 


---- 秦岭 边界 


50 100km 


107° 108。 109°E 
图 1 采样 点 及 气象 站 位 置 示意 图 


Fig.l Locations of tree ring sample’ site and meteorological 


stations 


主 序 列 相 关 性 差 的 个 别 树 芯 也 要 去 除 掉 ,最 终 剩 下 46 
条 轮 宽 序列 。 
1.2 太白 红 杉 树 轮 宽度 年 表 建 立 

由 于 研究 目的 为 重建 历史 气候 序列 ,所 以 在 完成 交叉 定年 后 ,要 剔除 轮 宽 序列 中 的 非 气 候 信 号 。 年 表 的 
建立 采用 ARSTAN 程序 "完成 ,在 此 过 程 中 ,采用 负 指 数 函 数 或 线性 函数 对 各 宽度 序列 进行 拟 合 ,最 后 得 到 
标准 化 年 表 ( STD) 和 差 值 年 表 ( RES) ,为 尽量 保存 年 表 的 低频 信息 ,采用 标准 化 年 表 进 行 气候 响应 和 重建 分 
析 。STD 年 表 的 统计 特征 列表 见 表 1 。 

年 表 的 起 始 年 代 以 子 样本 信和 号 强度 SSS( Sub 一 sample signal strength ) [8] 来 确定 ,为 尽 可 能 延长 年 表 , 本 研 
究 取 SSS 门限 为 0.75 ,结合 Rbar 值 , 确 定 可 售 年 表 时 段 为 1870 一 2013 年 。 从 表 1 可 以 看 出 ,STD 年 表 平 均 敏 
感度 (MS ) 较 高 ,反映 树 轮 宽窄 变化 大 人 径 网 生长 对 环境 变化 的 敏感 度 较 高 。 树 间 相 关系 数 、 信 噪 比 (SNR) 、 第 
一 主 成 分 方差 解释 量 (PC1) 及 样 杰 总体 代表 性 (EPS ) 均 较 高 ,表明 采样 点 树 轮 宽度 变化 有 较 好 的 一 致 性 。 


表 1 采样 点 概况 及 树 轮 样本 的 统计 特征 


Table 1\， Generalsituation of sampling sites and the statistical characteristics of tree ring samples 


特征 参数 Statistical charaateristiGs STD 特征 参数 Statistical characteristics STD 
采样 海拔 Elevation of_samplihg site8/m 3235 树 内 相关 系数 Correlation within trees 0.60 
所 用 样本 量 Samples used/( 棵 Y 芯 ) 26/46 “| 树 间 相关 系数 Correlation between trees 0.48 
年 表 平 均 敏 感度 Meansensitivity 0.19 样本 总 体 代表 性 Expressed population signal 0.92 
年 表 标准 差 Standard deviation 0.21 序列 长 度 Series length/a 165 
一 聊 自 相关 系数 )1*" okder auto correlation 0.51 子 样本 信号 强度 SSS>0.75( beginning year/ trees) 1870/9 
第 一 主 成 分 所 占 方 差 量 Variance in first eigen vector 54.04% | 子 样本 信号 强度 SSS>0.80 ( beginning year/trees) 1882/12 
信 品 此 , Signal to noise ratio 10.44 


1.3 ”气象 资料 处 理 

研究 区 由 于 海拔 过 高 ,附近 缺乏 长 期 的 器 测 数 据 ,为 实现 对 林 线 区 气候 变化 的 重建 ,我 们 利用 Arcgis 空间 
分 析 功能 通过 以 下 步骤 实现 多 元 数据 的 太白 山林 线 区 气温 资料 提取 : 

(1) 计 算 太 白山 逐 月 气温 直 减 率 。 根 据 太 白山 自然 保护 区 以 内 的 11 个 高 山 气象 站 点 的 2012 一 2014 年 气 
温 实测 资料 ,计算 太白 山南 , 北 坡 逐 月 的 气温 垂直 递减 率 , 结 果 见 表 2 . 

(2) 大 白山 空间 气温 搬 值 。 根 据 太白 山 自然 保护 区 以 外 的 太白 县 、 佛 坪 县 `. 眉 县 等 15 个 临近 气象 站 
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1983 一 2013 年 的 气温 资料 ,结合 来 自 国家 测绘 局 的 25 mx25 m 分 辨 率 DEM 高 程 数据 ,采用 普通 Kriging 插值 
法 !2 对 气温 要 素 进行 空间 插值 ,变异 函数 均 选 择 球面 函数 。 首 先 ,利用 表 2 中 的 垂直 递减 率 将 各 站 点 实测 气 
温 值 订正 至 海平 面 气温 ,在 ArcGIS10.1 中 进行 插值 计算 。 随 后 ,利用 栅 格 计算 工具 ,基于 气温 垂直 递减 率 将 
插值 后 的 栅 格 数据 与 区 域 DEM 数据 相 受 加 ,最 终 得 到 具有 地 形 特征 的 太白 山 保 护 区 整体 .保护 区 南北 坡 的 
气温 栅 格 数据 。 通 过 交 又 验证 显示 插值 结果 的 平均 标准 差 为 0.877, 对比 插 值 计 算 结果 与 部 分 高 山 气 象 站 点 
的 实测 数据 ,显示 插值 误差 在 +1°%C 以 内 ,表明 插值 结果 可 靠 性 较 好 ,可 以 代表 太白 山 保护 区 高 山 气候 特征 名 | 。 


表 2， 太 和 白山 高 海拔 区 南北 坡 逐 月 气温 垂直 递减 率 /(%/100 m) 


Table 2 ” Temperature lapse rate of Mt. Taibai alpine area from northern and southern aspects 


1 月 2 月 3 上 月 4 月 5 月 6 月 7 月 8 月 9 月 10 月 LE 朋 了 2 月 
January ”February March April May June July August September Qctober November WD JDecember 
北城 (CA100 m) 0.38 0.33 0.58 0.54 0.63 0.62 0.68 0.59 0.47 0.48 0i47 0.44 
North aspect 
南 坡 ( CA100 m) 0.44 0.55 0.56 0.56 0.50 0.57 0.54 0.55 0.55 0.48 0.55 0.48 


South aspect 


(3) 提 取 林 线 31 a 逐 月 平均 气温 。 采 用 点 值 提取 的 方法 ,在 具有 地 形 特征 的 太 自 山南 北 坡 气温 栅 格 数 据 
中 提取 了 海拔 为 3235 m 的 采样 区 (N 34°00'00,E107°42')1983 一 2013 年 的 逐年 逐 月 气温 ,并 计算 这 31 a 的 
1 一 12 月 的 逐 月 平均 气温 值 7 。 

将 步骤 (1) 中 距离 采样 点 距离 最 近 的 汤 峪 2 号 (海拔 :2767 my 经 纬度 ;N34°00'01,E107°48'56”) 高山 气 象 
站 点 的 气温 实测 数据 与 7 进行 比较 (图 2) ,可 见 , 逐 月 平均 气温 变化 趋势 吻合 的 同时 ,两 地 温差 也 符合 该 地 
区 的 气温 直 减 率 特征 。 考 虑 到 降水 随 海拔 分 布 特征 没有 明显 的 规律 , 且 受 地 形 影 响 较 明显 ,因此 插值 方法 暂 
未 应 用 于 降水 数据 ,而 是 采用 太白 山 周 边 的 太白 、 佛 坪 两 地 月 总 降水 量 的 平均 值 ,两 地 降水 数据 如 图 3 所 示 。 
可 见 , 太 白山 林 线 最 高 平均 气温 出 现在 7 月 为 12.05% ,而 最 低 为 12 月 的 -11.86C ,这 与 海拔 相对 较 低 的 太 
白 、 佛 坪 两 气象 站 月 均 温 的 最 高 值 (20.81% ) 及 最 低 值 ( -1.83% ) 存 在 明显 差异 "| 。 


1 下 200 『 
-@ 平均 气温 提取 值 Th 180 上 佛 坪 降水 
3 10 上 一 汤 峪 2 号 和 160 上 国王 太白 降水 
oO 
中 5 全 140 上 
本 号 120 上 
5 0 号 100 二 
号 12 £ 80 上 
路 “| 关 00[ 
也 -10 上 壮 40 上 1 
20 上 
-1$ 上 0 [mw ml | i | ml 1 
4 5 


1 2 3 6 - 8 9 ‘10 11 12 
图 2 采样 点 深 月 平均 气温 提取 值 ( 7, ) 与 汤 峪 2 号 气象 站 平均 气 月 份 Month 
人 图 3 太白 、 佛 坪 逐 月 平均 总 降水 量 分 布 


TEACogyisor 户 between extracted monthly mean air Fig.3 Monthly mean total precipitation of Taibai county and 


temperature ( 7，) and instrumental monthly mean air Foping county 


temperature from Tangyu No.2 meteorological station 


1.4 太白 红 杉 径 向 生长 主导 气候 因子 的 提取 

基于 上 文 处 理 的 气温 \ 降 水量 及 不 同月 份 组 合 ,考虑 到 气候 要 素 对 植物 生长 的 滞后 作用 ,选择 上 年 9 月 至 
当年 10 月 为 分 析 时 段 与 树 轮 宽 度 指数 序列 进行 Pearson 相关 分 析 。 

相关 分 析 结 果 显 示 ,太白 红 杉 与 气温 的 相关 性 较 与 降水 的 相关 性 强 且 大 致 呈正 相关 (图 4) 。 序 列 与 当 
年 生长 季 2 一 7 月 逐 月 气温 呈 不 同 程度 的 正 相 关 ,并 在 3.4.5,.6 月 均 通 过 了 95% 的 显著 性 检验 。 将 当年 生长 


http://www.ecologica.cn 


22 期 秦 进 


等 : 近 144 年 来 秦岭 太白 山林 线 区 3 一 6 


平均 气温 的 重建 7589 


显示 


季 1 一 9 月 不 同时 段 的 气温 组 合 后 ,相关 分 析 结 
序列 与 太白 山林 线 地 区 1 一 9 月 整体 平均 气温 的 相关 系 
数 达 0.519( P<0.001) ,尤其 与 3 一 6 月 平均 气温 的 相关 
值 达到 最 高 (r=0.671,P<0.001)。 而 序列 与 3 一 7 月 降 
水 均 呈 负 相 关 关 系 ,降水 负 相 关 最 强 的 时 段 与 气温 正 相 
关 最 强 的 时 段 相对 应 , 均 为 3 一 6 月 ,与 4 月 单 月 总 降水 
量 的 负 相 关 最 高 (r=0.398 ,P<0.05)。 反 映 研究 区 生长 
季 气 温 一 定 程度 的 增高 是 有 助 于 树木 生长 的 ,有 利于 较 
宽 年 轮 的 形成 ,但 生长 季 降 水 若 增多 则 会 对 树木 生长 产 
生 明 显 抑制 作用 。 因 此 ,3 一 6 月 气温 是 太白 红 杉 年 轮 
生长 的 主导 因子 。 

从 树木 生理 学 角度 来 看 ,树木 年 轮 的 宽度 波动 往往 
是 由 一 个 或 多 个 气候 因素 共同 作用 的 结果 ,而 本 研究 样 
本 均 采 自 3000 m 海拔 以 上 的 极端 生境 ,树木 进行 正常 
生命 活动 所 依赖 的 气温 与 降水 都 较 荐 乏 ,但 是 ,秦岭 地 
区 具有 海拔 高 .气候 湿润 . 土 层 较 厚 及 土壤 含水 量 高 的 
特点 '” ,在 树木 生长 季 前 期 ,气温 过 低 会 无 法 满足 树木 
维持 正常 生长 代谢 活动 的 需要 ,使 树木 形成 层 活 动 推 
迟 ,从 而 形成 较 守 的 年 轮 ,而 气温 升 高 有 利于 地 温 升 高 ， 
促使 植物 根系 开始 活动 , 当 气温 升 高 至 树木 生长 的 最 低 


时 气温 
号 降水 


相对 系数 Correlation coefficients 


9 1011 12 1 2 3 4 5 6 7 8 9 103 一 0 
月 份 Month 


图 4 树 轮 宽 度 STD 年 表 与 逐 月 平均 气温 及 总 降水 量 相关 结果 
Fig.4 The correlation analysis between STD chronology and 


monthly mean temperature ,total precipitation 
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温度 时 ,树木 形成 层 细胞 开始 开裂 ,有 利于 延长 生长 季 
和 星 材 宽 年 轮 的 形成 。 基 于 以 上 原因 ,在 秦岭 主 
疹 太 白山 林 线 地 区 ,春季 至 初夏 的 气温 成 为 树木 生长 的 
主要 限制 因子 ,而 前 人 在 研究 秦岭 其 他 地 区 树种 


图 5 3 一 6 月 气温 实测 值 与 重建 值 的 去 趋势 化 指数 比较 
Fig.S The comparison between the detrended indices of 
reconstructed temperature ( March 一 June ) and the instrumental 


data 


的 过 程 中 也 得 到 了 生长 季 气 温 为 树 本 生长 主导 因子 的 结论 。 


2 结果 与 分 析 


2.1 回归 模型 的 建立 与 调整 


根据 相关 性 计算 结果 和 以 太白 山林 线 地 带 3 一 6 月 气温 作为 因 变 量 , 太 白山 太白 红 杉 STD 年 表 作 为 自 变 
量 ,采用 线性 回归 方法 对 太白 山林 线 公 元 1870 至 2013 年 以 来 3 一 6 月 平均 气温 进行 模拟 重建 ,转换 方程 为 : 


7T,, = 2.505ISTD + 3.932 


(1) 


式 中 ,7 名 为 3 守 6 和 绷 平 均 气 温 ,1sy 为 太白 红 杉 STD 年 表 , 该 方程 的 相关 系数 r=0.670(P<0.001) ,方差 解释 量 
为 44.9% , (调整 生 由 度 后 方差 解释 量 为 42.8%) ,=21.184。 鉴 于 从 插值 所 得 气象 记录 只 有 31 年 ,我 们 采用 
逐一 浏 | 除法 来 检验 转换 方程 的 稳定 性 。 发 现存 在 奇异 年 份 1997 .1998 年 使 方程 (1) 稳定 性 下 降 , 去 除 1997、 
1998 年 后 ,相关 系数 + 升 为 0.757(P<0.001)。 气 象 观 测 记录 表明 ,1996 年 出 现 了 3 一 6 月 气温 均 低 于 多 年 平 
均值 的 现象 ,其 滞后 作用 可 能 对 1997 树木 年 轮 组 织 生长 造成 了 较 显 著 的 影响 ,虽然 1998 年 3 一 6 月 气温 有 所 
回升 ,但 是 3 一 5 月 降水 却 出 现 了 多 年 来 的 罕见 极 高 值 (84.95 mm) ,树木 生长 受到 严重 抑制 ,产生 罕 轮 ,加 之 
1997 1998 年 我 国 气候 受 ENSO 事件 严重 波及 发 生 异 常 .中 ,共同 导致 回归 模型 对 1997 .1998 年 3 一 6 月 平均 
气温 的 拟 合 结果 偏差 较 大 ,因此 考虑 从 回归 方程 中 剔除 '” 。 在 去 除 1997 .1998 年 之 后 ,转换 方程 为 : 
Ts = 3.015ISTD + 3.355 (2) 
由 表 3 可见, 与 方程 (1) 相 比 ,方程 (2) 在 各 个 统计 指标 上 优势 更 明显 ,表明 方程 更 加 稳定 。 经 过 交叉 检 


http://www.ecologica.cn 


7590 生 态 学 报 37 卷 


验 分 析 ,重建 序列 与 观测 序列 的 符号 检验 (S1) .一 阶 差 符号 检验 (S2 ) ,误差 缩减 值 ( RE ) ,乘积 平均 数 (1) 以 及 
方程 的 下 检验 值 等 统计 量 均 体 现 了 重建 结果 的 可 靠 性 。 图 5 为 运用 转换 方程 (2) 重建 的 太白 山林 线 3 一 6 月 
平均 气温 与 实测 值 插值 结果 的 对 比 。 在 观测 时 段 内 (1983 一 2013 年 ) ,去 趋势 后 的 重建 值 与 实测 值 吻合 较 好 ， 
证 明 重建 结果 可 靠 , 能 够 反映 太白 山林 线 地 区 气温 波动 规律 '”| 。 


表 3 各 项 检验 参数 


Table 3 Calibration and verification statistics 


时 段 Period r R? Ra Sl S2 RE F t p 
1983 一 2013 0.757 57.2% 55.5% 20(19* ,21** ) 23(21* ,22**) 0.543 32.129 2.571 <0.001 
r 为 相关 系数 , 尼 为 方差 解释 量 ,Ra 为 调整 自由 度 后 的 方差 解释 量 ,S1 为 符号 检验 ,S2 为 一 阶 差 符 号 检验 , RE 为 误差 缩减 值 ,i 为 乘积 平均 
数 , 为 方差 检验 值 ; * 代表 95% 的 置信 水 平 , * * 代表 99% 的 置信 水 平 


2.2 太白 山林 线 区 3 一 6 月 平均 气温 的 重建 
基于 转换 方程 (2) ,将 子 样本 代表 性 >0.75 的 年 份 1870 年 作为 起 始 年 ,2013 作为 涤 止 窜 玫 R 们 重建 了 太白 
”山林 线 地 区 历史 时 期 3 一 6 月 平均 气温 (图 6) 。 
与 可 见 ,1870 年 以 来 ,研究 区 3 一 6 月 平均 气温 波动 频繁 ,波动 幅度 在 1.83 一 4.85 也 过 间 ,最 低 和 最 高 温度 相 
差 达到 3.06%C ,重建 时 段 3 一 6 月 气温 平均 值 7 ,为 3.51% ,标准 差 合 为 463% 在 本 征 究 中 ,定义 气温 高 于 
7T ,+1o(4.14%C ) 的 年 份 为 高 温 年 份 ,气温 低 于 7, -10(2.88%C ) 的 年 份 为 低温 年 份 , 而 介 于 两 者 之 间 的 气温 
”为 正常 值 ”。 由 表 4 可 见 ,在 过 去 144 年 中 , 偏 暖 年 份 出 现 19 次 偏 次 年 份 出 现 26 次 ,分 别 占 总 年 份 的 13% 
和 18%。 


[ 


~ > 一 一 重建 从 ”一 一 1 年 滑动 平均 信 


Theant+ GO 


| | 人 
;| NW 7 UN rv 2 i 
(i 


气温 Air temperature/°C 


样本 量 Sample depth 


1870 1890 1910 1930 1950 1970 1990 2010 
年 份 Year 


图 6 太白 山林 线 区 1870 一 2013 年 3 一 6 月 平均 气温 重建 序列 
Fig.60” The March 一 June temperature reconstruction for the timberline of Mt. Taibai from 1870 to 2013 


重建 结果 中 出 现 的 偏 暧 、 偏 冷 年 份 ( 表 4) 与 陕西 地 情 资料 库 近 200 年 所 发 生 的 气象 灾害 记录 有 和 较 好 的 
对 应 (数据 来 源 于 陕西 地 情 网 http ://www.sxsdq.cn ) 。 偏 暖 年 份 基本 伴随 着 大 范围 的 严重 干旱 事件 并 对 农业 
产生 了 严重 影响 ,例如 1892 .1928 .1929 .1959 年 ,记载 有 春 旱 . 夏 早 ,大 旱 等 现象 ,这 些 年 份 “收成 歉 薄 ”““ 粮食 
缺乏 ”““ 夏 秋收 成 不 到 二 成 ”。 而 偏 冷 年 份 在 陕西 出 现 了 较 大 范围 的 洪水 灾害 ,例如 1881 、1921 、1948 、1949 
年 ,记载 有 “ 雨 蜀 伤 禾 ”山水 陡 涨 ”河水 暴涨 "等 现象 其 中 连续 超过 2 a 及 以 上 的 高 温 时 段 为 1966 一 1968 
年 .2005 一 2007 年 ,连续 超过 2 a 及 以 上 的 低温 时 段 为 1871 一 1873 年 1880 一 1883 年 .1997 一 1999 年 。 

为 了 反映 重建 序列 的 低频 特征 ,更 好 地 观察 3 一 6 月 平均 气温 的 年 代 际 变化 ,对 序列 作 11 a 滑动 平均 ,并 
定义 滑动 平均 值 连续 11 a 以 上 高 于 平均 值 的 时 间 段 为 暖 期 ,滑动 平均 值 连续 11 a 以 上 低 于 平均 值 的 时 间 段 
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为 冷 期 2 。 由 图 6 可 见 , 在 重建 期 间 气温 偏 高 的 (>3.51%C ) 时 段 主要 有 1882 一 1892 年 .1919 一 1929 年 和 
1997 一 2013 年 ,总 持续 时 间 为 29 年 ,而 在 重建 期 间 气温 较 低 (<3.51% ) 的 时 段 主要 :1870 一 1881 年 、1903 一 
1918 年 和 1977 一 1996 年 ,总 持续 时 间 为 39 年 。 重 建 时 段 的 低频 变化 趋势 中 偏 冷 时 有 段 较 偏 暧 时 段 多 ,这 与 极 
端 年 所 占 总 年 的 比例 一 致 ,其 中 1977 一 1996 年 冷 期 持续 时 间 最 长 (22 a) ,20 世纪 90 年 代 以 前 咀 期 持续 时 间 
较 短 升温 幅度 较 小 。 总 体 而 言 ,19 世纪 末期 是 一 个 以 1881 年 为 界限 的 显著 冷暖 交替 时 期 ,在 20 世纪 初 至 
20 年 代为 一 显著 冷 期 ,而 20 一 30 年 代 末 是 一 个 显著 的 暧 期 ,1919 一 1931 年 的 高 温 期 与 我 国 北方 20 年 代 出 现 
的 干 量 事 件 相 吻合 '“ 。 随 后 ,在 1931 一 1978 年 这 一 时 期 ,40 a 的 初春 温度 变化 不 存在 明显 的 变 暖 趋势 号 温 
变化 幅度 较 小 ,冷暖 交替 频繁 ,气温 相对 比较 稳定 ,而 在 20 世纪 70 年 代 末期 气温 开始 下 降 ,从 20 世纪 90 托 
代 后 期 至 今 气温 又 明显 持续 增 暖 。 


表 4 144 年 以 来 偏 暖 、 偏 冷 年 份 及 重建 温度 值 


Table 4 Extreme years and reconstructed temperature in recent 144 years 


序号 。 偏 暖 年份 。 重建 值 C 售 淮 年 份 、 覃 建 值 /AC | 序号。 。 偏 暖 年份。 一 建 但 人 区 y 年 份 ”重建 值 / 
No1 We Reconstructed Golllyea Reconstructed Nol Wik yor Recongstructed CR year Reconstructed 
value value value value 
1 1886 4.636125 1871 2.679565 14 1977 4.490705 1949 2.51762 
2 1891 4.398165 1872 2.507705 15 1995 4.17873 1956 2.81507 
3 1892 4.93688 1873 2.758885 16 2005 4.36842 1964 2.87456 
4 1928 4.57333 1880 2.415165 17 2006 4.79146 1969 2.884475 
5 1929 4.89061 1881 1.826875 18 2007 4.26266 1971 2.765495 
6 1935 4.27588 1882 2.540755 19 2013 4.391555 1986 2.706005 
2 1943 4.206475 1883 2.692785 20 1987 2.580415 
8 1945 4.669175 1899 2.79524 21 1989 2.421775 
9 1959 4.77163 1920 2.216865 22 1990 2.249915 
10 1960 4.629515 1921 1:998735 23 1997 2.52423 
11 1966 4.24283 1933 2.765495 24 1998 2.276355 
12 1967 4.199865 1940 2.86795 25 1999 2.34576 
13 1968 4.325455 1948 2.29949 26 2010 2.818375 


2.3 重建 结果 的 对 比 与 验证 
2.3.1 重建 结果 冷暖 期 的 区 域 对 比 

通过 与 田 沁 花 等 '" ,所 重建 的 过 去 134 a 伏 牛 山 5 一 7 月 平均 最 高 气温 变化 趋势 进行 对 比 ,发现 该 研究 
中 得 到 的 偏 冷 时 段 1900 一 1917 年 、1974 一 1989 年 和 偏 暖 时 段 1889 一 1899 年 、1918 一 1940 年 、1990 一 2002 年 
与 本 研究 重建 所 得 偏 冷 期 与 偏 暖 期 基本 重合 ,然而 ,在 1931 一 1978 年 这 一 时 期 40 a 伏 牛 山 5 一 7 月 平均 最 高 
气温 变化 存在 着 明显 的 变 暖 趋势 ,但 在 相同 时 段 ,太白 山林 线 区 3 一 6 月 平均 气温 上 升 趋势 不 明显 。 与 刘 珊 、 
刘 娜 等 重建 秦岭 中 段 分 水 岭 地 区 1 一 7 月 平均 气温 的 结果 进行 比较 发 现 ,太白 山 与 佛 坪 光 头 山 均 在 在 19 世 
纪 末 纳 及 20 蛋 纪 30 年 代 之 前 存在 低 一 高 一 低 一 高 的 两 次 显著 震荡 ,在 1934 一 1990 年 这 一 时 期 的 气温 变化 
均 不 存在 明显 的 变 暧 趋势 ,气温 变化 幅度 小 ,气温 相对 稳定 ,持续 到 20 世纪 90 年 代 初 期 又 出 现 了 明显 的 气温 
升 高 现象 。 而 蔡 秋 芳 .刘备 '31 等 ,在 陕西 中 北部 地 区 1826 年 以 来 4 一 9 月 气温 的 重建 结果 中 发 现 ,20 世纪 20 
年 代 末 期 是 自 1826 年 以 来 陕西 中 北部 地 区 4 一 9 月 温度 最 高 的 时 段 ,陕西 中 北部 地 区 自 1923 年 以 来 4 一 9 月 
平均 气温 不 断 上 升 ,到 1929 年 温度 达到 重建 时 段 内 最 高 记录 ， 而 在 本 研究 区 ,20 世纪 20 年 代 末 期 同样 为 
3 一 6 月 气温 的 最 高 时 段 ,并 且 , 重 建 气温 自 进入 20 世纪 以 来 ,极端 最 高 气温 也 出 现在 1929 年 (4.89% ) ,而 这 
次 极端 高 温 事件 所 导致 的 灾害 在 中 国 北方 其 他 地 区 的 树 轮 研究 中 也 得 到 了 印证 5 5 。 

与 甘肃 省 东南 部 天 水 石门 山 5 一 6 月 最 高 气温 的 树 轮 重建 结果 (数据 来 源 于 http://www. ncede. noaa. gov/ 
paleo/study/17033 ) 进行 比较 (图 7) 发 现 ,两 研究 区 高 、 低 温 波动 大 体 一 致 ,相关 系数 为 0.171 (N= 133,P< 
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0.05) ,表明 在 低频 变化 上 两 序列 是 显著 相关 的 。 冷 暖 时 
段 的 变化 上 两 者 出 现 局 部 时 段 的 不 同步 性 ,这 可 能 是 由 于 
所 表征 的 月 份 、 气象 要 素 以 及 海拔 并 不 完全 一 致 造成 的 ， 
然而 20 世纪 30 年 代 的 高 温 期 ,40、50 年 代 的 低温 期 .60、 
80 年 代 的 高 温 期 以 及 90 年 代 以 来 的 气温 明显 上 升 趋势 在 
两 条 曲线 中 共存 ,表明 太白 山林 线 地 带 与 石门 山地 区 的 气 
温 变 化 序列 经 历 了 极为 相似 的 冷暖 波动 过 程 ,两 个 曲线 的 
同步 变化 表现 出 较 好 的 区 域 一 致 性 ,印证 了 本 重建 结果 的 
可 靠 性 。 
2.3.2 太白 山林 线 区 3 一 6 月 平均 气温 的 周期 性 

本 研究 选用 Morlet 连续 小 波 变换 对 气温 序列 进行 分 
析 ,为 减弱 小 波 变 换 资料 的 始 . 未 区 域 所 对 应 变换 结果 的 
边界 效应 ,文中 采用 对 称 延伸 法 减弱 这 种 影响 '*i。 结 合 图 
6 可 清楚 看 到 太白 山林 线 地 区 3 一 6 月 平均 气温 在 过 去 144 
a 的 演化 过 程 中 存在 多 时 间 尺 度 特征 并 与 方差 图 中 的 峰 谷 
变化 一 致 。 与 重建 序列 曲线 (图 7) 进行 比较 可 以 发 现 , 周 
期 信号 比较 强 的 时 段 气 温 波 动 也 相对 明显 ,周期 信号 不 明 
显 的 时 段 ,气温 的 波动 相对 缓和 。 

太白 山林 线 区 3 一 6 月 气温 演变 过 程 中 存在 着 相对 显 
著 的 31 一 22 a,22 一 18 a 以 及 13 一 10 a 的 3 类 尺度 准 周 期 


27.5 
27.0 
26.5 
26.0 
25.5 
25.0 
24.5 


24.0 — 
5.0 


4.5 
4.0 
3 
3.0 
2.5 
2.0 


气温 Air temperature/°*C 


1.5 
1870 1890 1910 1930 1950 1970 1990 2010 
年 份 Year 


7 同期 石门 山 5 一 6 月 最 高 气温 (a) 与 太白 山林 线 3 一 6 月 
平均 气温 (bjt 

Fig.7A Comparisor between the reconstruction of May 一 June 
maximum temperature of Mt. Shimen (a) and the 
reconstruction of March 一 June timberline mean temperature 


of Mt. Taibai (b) at their interval 


变化 规律 。13 一 10 a 尺度 的 周期 信号 在 时 间 尺 度 .上 并 不 稳定 9 ,振幅 随时 间 推 移 产 生 了 明显 变化 的 情况 ,在 
1900s 之 前 和 1940s 以 后 表现 的 较为 稳定 ,但 礁 1900s 一 940s 和 1990s 一 2010s 期 间 发 生 了 较 大 浮动 。22 一 18 
a 的 变化 周期 与 马 利 民 、 刘 盘 '” 等 ,对 秦岭 SSOI(The Southern Oscillation Index) 重 建 序列 的 功率 谱 分 析 中 发 现 
的 自 1743 年 以 来 该 序列 存在 明显 的 21. 82 一 17.14 a 周期 比较 吻合 ,而 22 年 的 周期 是 Hale 周期 ,这 很 可 能 说 
明了 南方 涛 动 对 太白 山林 线 区 树林 生长 存在 显著 影响 。 在 13 一 10 a 尺度 上 出 现 的 显著 冷 一 暖 交 替 准 17 次 
震荡 ,与 太阳 黑子 活动 周期 特征 相似 并 对 应 了 西北 干旱 灾害 的 11 a 准 周 期 '*。 这 些 周期 暗示 全 球 的 气候 变 
化 特征 也 存在 于 该 地 区 的 气候 房 虹 ,太白 山林 线 地 区 3 一 6 月 气温 可 能 受到 了 大 范围 的 气候 变化 影响 。 
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图 8 太白 山林 线 区 气温 重建 序列 的 小 波 分 析 


Fig.8 Wavelet analysis for the reconstructed temperature series of Mt. Taibai timberline 


http://www.ecologica.cn 


201712.01484v1 


图 
国 


chinaXiIv 


ChinaXiv 合 作 期 刊 


22 期 秦 进 等: 近 144 年 来 秦岭 太白 山林 线 区 3 一 6 月 平均 气温 的 重建 7593 


3 结论 


论文 以 太白 山林 线 区 的 太白 红 杉 树 轮 样 本 建立 的 STD 年 表 为 基础 ,与 借助 Arcgis 软件 气温 插值 功能 提 
取 的 采样 点 过 去 31 a 气温 数据 进行 相关 分 析 , 发 现 树 轮 年 表 对 温度 变化 响应 敏感 ,与 研究 区 生长 季 气 温 呈 
著 正 相关 。 在 此 基础 上 ,实现 了 太白 山林 线 1870 一 2013 年 间 3 一 6 月 平均 气温 的 重建 ,并 得 到 如 下 结论 : 

(1) 重 建 结 果 显 示 ,太白 山林 线 区 3 一 6 月 平均 气温 在 过 去 144 a 经 历 了 3 次 冷 期 和 3 次 暧 期 的 波动 19 
世纪 末 与 20 世纪 初 的 冷 一 暧 交替 ,20 世纪 20 一 30 年 代 末 的 显著 高 温 期 1931 一 1978 年 间 相 对 比较 稳 宠 的 冷 
暖 波 动 以 及 在 20 世纪 90 年 代 初期 出 现 的 明显 气温 升 高 等 现象 与 周边 地 区 同期 气温 重建 结果 相似 。 

(2) 研 究 区 3 一 6 月 平均 气温 重建 序列 与 天 水 石门 山地 区 5 一 6 月 最 高 气温 同期 重建 结果 在 高 低频 变化 
上 均 表 现 出 较 好 的 一 致 性 ,从 另 一 方面 佐证 了 本 研究 方法 以 及 重建 结果 的 可 靠 性 。 

(3) 周 期 分 析 结 果 表 明 ,太白 山林 线 地 带 气温 存在 31 一 22 a,22 一 18 a 以 及 13; -108 的 引种 主要 变化 周 
期 ,周期 特征 与 太阳 活动 及 SOI 的 周期 特征 接近 ,意味 着 太白 红 杉 宽度 可 能 包含 更 大 局 度 的 气候 波动 信号 。 
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验 提供 了 LINTAB 树 轮 分 析 仪 和 实验 室 , 以 及 包 光 老师 在 对 年 轮 梯 品 的 处 理 及 分 析 过 程 中 给 予 的 极 大 技术 
帮助 。 
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